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We review the current theory of atom lasers. A tutorial treatment of second 
quantisation and the Gross-Pitaevskii equation is presented, and basic concepts 
of coherence are outlined. The generic types of atom laser models are surveyed 
and illustrated by specific examples. We conclude with detailed treatments of the 
mechanisms of gain and output coupling. 

1 Introduction 

An atom laser is a device which produces a bright, directed, and coherent 
beam of atoms. An ideal atom laser beam is a single frequency de Broglie 
matter wave, approximating a noiseless sinusoidal wave; in particular it will 
have a well defined frequency, phase, and amplitude. Quantum mechanically, 
such a wave is a mode of the quantised field, and the properties of brightness 
and well defined phase require the mode to be highly occupied. This is equiv- 
alent to reauiring Bose degeneracy, and rules out the possibility of a fermionic 
atom laserEI. 

Atom lasers are a recent concept, and to our knowledge the earliest scien- 
tific description appeared in 1993 u. ^Xhs first comprehensive papers on matter 
wave amplificationu and atom lasersQQQ appeared-in 1995, and the first exper- 
imental demonstration of an atom laser in 1997 u. The essential components 
of an atom laser have been identified by analogy with optical lasers and are 
illustrated in Fig. 1. They are; a reservoir of atoms from which the laser is 
pumped, a laser mode in which atoms are accumulated, a pumping (or gain) 
mechanism which transfers atoms from the reservoir into the laser mode by 
stimulated transitions, and an output coupler which produces an output beam 
from the laser mode while preserving its coherence. 

Many of the possible applications of atom lasers are probably still to be 
imagined, nevertheless atom lasers have the potential to revolutionise atom 
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Figure 1. Schematic diagram of a generic atom laser. 



optics just as optical lasers did conventional opticaS'Ll. A bright coherent beam 
greatly enhances interferometric capabilities, for example allowing fringes to 
be obtained on a small time scale, and the use of unequal path length in- 
terferometers. The wavelength of the atoms is typically much smaller than 
light, and an obvious potential technological application is the precision de- 
position of matter for nanofabrication. While the similarities to optical lasers 
are important, the fundamental differences between atoms and photons are 
also significant. Atoms interact with each other, are subject to gravity, have 
complex internal structure, and may travel at any speed, not just c. Their in- 
teractions (collisions) make atom optics intrinsically nonlinear, as illustrated 
by a recent experiment showing self-induced four-wave mixing of atom laser 
beamalj. Their interaction with the gravitational field gives the opportunity 
for atom laser beams to be used in high precision gravitational and iner- 
tial measurementsllj, and verv sensitive atom interferometers may be able to 
detect changes in space-timeQ. Eventually we expect that atom lasers will 
provide a convenient source for greatly improved atomic clocks, or for manip- 
ulation of the quantum mechanical state of the matter wave, analogous for 
example to optical squeezingla. 

As of April 2000 five experimental grouijisJiajte demonstrated atom lasers 
experimentally, in either pulsed operationEHlil'tj or with quasi-continuous 
outputll3'tZl. In each case, the laser mode consists of a Bose-Einstein con- 
densate confined in the ground state of a magnetic trap, gain is achieved by 
evaporatively cooling a reservoir of thermal atoms, and the output coupling 
is accomplished by using an external electromagnetic field to transfer the in- 
ternal state of the atom from a trapped to an untrapped state. In these 
experiments the output coupling and the pumping processes are separated in 
time. Q-switching in pulsed optical lasers also separates the pumping and 
output coupling, but the cycle time is much shorter. Hence, a closer anal- 
ogy is to a cavity loaded with coherent light which then leaks out through a 
mirror. The use of the name "atom laser" is justified by the coherence of the 
output beam, which was first demonstrated by the MIT groupaE3, and is due 
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fundamentally to the stimulated emission of bosonie atoms into the ground 
state of the trap. It remains an outstanding experimental goal to achieve 
simultaneous pumping and output coupling to make a truly continuous atom 
laser. 

Despite the parallels with optical laser theory, the fundamental differences 
between photons and atoms remain important. They arise principally because 
atoms have rest mass, and because they interact with each other. Unlike 
photons, atoms cannot be created or destroyed, rather they are transferred 
into the laser mode from some other mode. The interactions produce phase 
dynamics, which degrades the coherence of the atom laser, and self repulsion 
which spreads the output beam and limits the possible focussing. Even when 
these interactions are neglected, atoms have dispersive propagation in vacuum. 

Various critical parameters for atom and optical lasers arc radically dif- 
ferent. For example the characteristic frequency in the optical case is the 
frequency of the light, of order 10^^ Hz, while for atoms the de Broglie fre- 
quency is closer to that of the trapping potential, of order lO'' Hz. And, 
whereas it is so far only practical to get atoms to accumulate into the ground 
mode of a trap, optical lasers typically operate on high modes of the optical 
cavity. Finally, a practical consideration is that atom lasers can operate only 
in ultra high vacuum. 

In this paper we review the current state of theory of atom lasers. We 
begin, in sections 2-4, with a brief tutorial review of second quantisation and 
the Gross-Pitaevskii equation, and the concepts of coherence. In section 5 
we survey the main types of theoretical models for atom lasers, and illustrate 
these with specific examples. Finally, in sections 6 and 7 we provide more 
detailed treatments of the mechanisms of gain and output coupling. 

2 Second Quantization 

In this section we review the theory of quantum fields: the most fundamen- 
tal theorctic;al framework for treating atom lasers. As long as their internal 
atomic structure is not disrupted, atoms may be approximated as field quanta. 
However the internal state of the atom may change, e.g. between trapped and 
untrapped orientational spin states, and then each relevant internal state is 
described by its own field. 

Identity of particles arises naturally in quantum field theory, as all par- 
ticles of the same type are indistinguishable quantised excitations of the 
same field. We introduce quantum fields by second quantization of the 
wavefunction. For a more complete treatment see Quantum Mechanics by 
E. Merzbacher. The Hamiltonian for N particles interacting via the pairwise 
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potential V (r, r') is 

N 

" Hn (rA + ^ 



^ 1 

i/ = ^i/„(r,) + -5]y(r„r,), (1) 



where the are the particle position vectors. The factor of a half in the 
2-body interaction term corrects for double counting of particle pairs. The 
single particle Hamiltonian Hq (r) typically represents the kinetic energy and 
an external potential Ve (r), such as a trap. In the position representation it 
is 

Ho = -—\7^+VEir). (2) 
2m 

Second quantization allows us to automatically account for the identity of the 
particles, and is achieved by introducing the operator field \1/ (r) which anni- 
hilates a particle at position r. It obeys the following canonical commutation 
rules with its hermitian conjugate particle creation field (r), 

(r) , (r')] = (5 (r - r') , [* (r) , ^' (r')] = 0. (3) 

The second quantized Hamiltonian has the form of an "expectation value" of 
the Hamiltonian (|l|) with respect to these fields; 

H ^ j (r) Ho^ (r) d'V +]^ j j -^"^ {v) (r') V (r, r') * (r') ^' (r) dVdV. 

(4) 

In the Heisenberg picture, operators obey the usual Heisenberg equation of 
motion. For example the field operator obeys 

in^^^^[^'{v,t),H]. (5) 
Using the commutation rules (|^) and the Hamiltonian (^) this becomes 

= |i/o (r) + y (r') V (r, r') * (r') * (r) . (6) 

It is often useful to expand the field operators in terms of products of operators 
tti and the elements of a basis set of orthonormal mode functions (r)}, 

* (r) = ^ a,0, (r) ^ a, ^ j 4>* (r) ^ (r) d\ (7) 

This expansion is the basis of the Fock, or particle number, representation 
of the field. The commutation relations for the annihilation and creation 
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operators follow from Eqs.(^) and the mode function orthonormality, 

a„ 4] = J J^r (r) 0; (r') (r) , *t (r')] dVdV = <5,, (8) 

These commutators are those for a set of harmonic oscillators. Hence we can 
identify d'i 0jS cL harmonic oscillator type annihilation operator for the mode 
(pi, and its Hermitian conjugate a| as the corresponding creation operator. In 
quantum field theory these are interpreted as annihilating and creating field 
quanta, which in our case are atoms, in the mode. The quantum mechanical 
state of the system is that of these oscillators. In the Fock representation the 
second quantised Hamiltonian (Q) becomes 

H = J2 / ^* Ho (r) (r) (9) 
i J2 ala]akai f f <f>* (r) <f>* (r') V (r, r') 0^ (r') (r) d\d\' 



2 . . 

If the mode functions are chosen to be eigenstates of the single particle free 
Hamiltonian, i/o^j — Eicjji, the free part of the Hamiltonian can be expressed 
in terms of the oscillator number operators {a|^aj}. Denoting the integrals in 
the interaction part of the Hamiltonian by Vijki we then find 

H ^'^aluiEi + - ^ aja'jakaiVijki- (10) 

i ij,k,l 

Instead of solving for the full dynamics one may find the quasi-particles: the 
weakly interacting excitations of the system. The aim is to find an (approx- 
imately) diagonal form of H. The field is expressed in terms of the quasi- 
particle creation foj and annihilation bi operators by 

* (r) = ^ (u, (r) 6,; ~ V* (r) bf) + V-c (r) . (11) 

i 

ipc is called the condensate, or macroscopic, wavefunction. It is normalised 
so that the integral over all space of its squared modulus is Nc , the number 
of particles in the condensate. Assuming a non-zero condensale. wavefunction 
is the basis of the "spontaneous symmetry breaking" methodca for analysing 
BECs, which is discussed further in the next section. The functions Ui (r) and 
Vi (r) are determined such that the Hamiltonian ( p^ has the approximate and 
noninteracting form 

H ^l^^^'>. + constant (12) 
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The quanta of these modes are the quasi-particles. The transformation to a 
quasi-particle representation is known as the BogoUubov transformation. 

3 The Gross-Pitaevskii Equation 

In an atom laser the atoms will interact, both in the laser itself and in the 
output atomic beam. The interactions between cold atoms are usefully de- 
scribed as collisions. Elastic collisions conserve the total external energy of 
the colliding atoms. Inelastic collisions do not, and energy may be exchanged 
between the external and internal degrees of freedom. Changes in internal 
state due to inelastic collisions are an important source of losses of atoms 
from magnetic traps. We consider a simple s-wave scattering model of elastic 
collisions. This approximation is valid for low temperatures and low densi- 
ties, such that both the atomic de Broglie wavelength and the inter-particle 
separation exceeds the characteristic range of the interaction potential. For 
further discussion see the article by Burnett in this volume. 

In this limit the atom-atom potential may be approximated by a delta 
function potential with the same interaction energy 

V{r,r')^Uo6{r-r'), Uo ^ -. (13) 

m 

a is called the "scattering length" . It is positive for repulsive interactions and 
negative for attractive ones. For indistinguishable particles the scattering 
cross section (the ratio of scattered to incident particle flux) is cr = Sttc^, 
rather than 47ra^, as found for distinguishable particles. Uq is the atom-atom 
interaction energy per atom. Using this approximation in the field operator 
equation (|[) gives 

i^^^ = {Ho (r) + Uo^^ (r) (r)} <f (r) . (14) 

We next approximate this operator equation by a non-linear Schrodinger equa- 
tion for the condensate wavefunction. 

The condensate wavefunction ipc is defined by Eq.(^T]) to be the con- 
tribution to the expectation value of the field operator ^E" due to the quasi- 
particle vacuum state. According to the "spontaneous symmetry breaking" 
assumption, above the BRC transition temperature ipc = 0, while below this 
temperature it is nonzeroEa. The name arises because the Hamiltonian (m 
is symmetric, or invariant, under a global phase change 5* (r) ^ 4' (r) e , 
whereas the condensate wavefunction must have a particular phase, breaking 
the global phase symmetry. This symmetry is a consequence of the pairing 
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of particle creation and annihilation operators, and hence of particle num- 
ber conservation, in the Hamilton i an . (^). Particle number conservation in 
fact implies that ^(r) = alwayslialij, but spontaneous symmetry breaking 
is a useful heuristic assumption, which used with due caution, gives experi- 
mentally verified results. In fact, a number of calculations have been made 
without this assumption, and obtained the same results. Grifhn's article in 
this volume further explores the spontaneous symmetry breaking assumption. 

At zero temperature there are no quasi-particle excitations and the field 
operator ^ in Eq. (|l^) may be approximated by the single particle condensate 
wavefunction tpc- We further assume that products of field operators may 
be approximated by the corresponding products of condensate wavefunctipns. 
The condensate wavefunction then obeys the Gross-Pitaevskii equationEl 

^^d^^cjrl ^ / |lv2 + Ve (r) + Uo Wl'l (r) . (15) 

Mathematically, this equation is a non-linear Schrodinger equation about 
which nmch is known. Nevertheless, in general a numerical solution is 
requiredEil. 

4 Coherence 

An atom laser is a device which emits a coherent beam of atoms. Coherence 
and noise are closely related. A noiseless classical wave, that is a perfect 
single frequency sinusoid, is defined to be perfectly coherent or perfectly cor- 
related. For such a wave the field amplitude at any space-time event A{x, t) 
completely determines the amplitude at any other event. In practice the field 
amplitude fluctuates due to classical or quantum mechanical noise, and such 
determination is not possible. Classical noise may be technical, such as vi- 
brating, components, or dynamical, such as relaxation oscillations in optical 
laser£3. Quantum mechanical noise may be interpreted as a consequence of 
the uncertainty relations. 

Optical lasers operating far above threshold have a well stabilised inten- 
sity. Semiclassical theories describe optical lasers by quantising the atoms 
only, not the light. Introducing a phenomenological noise term allows sponta- 
neous emission to be modelled without quantising the light. For spch a theory, 
using the dimensionless intensity I defined by Mandel and Wolfed, far above 
threshold 

((A/)2)«2, ((A/)2)V2/(j)^(V2)/(/) (16) 
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The angle brackets indicate the average over the noise which models sponta- 
neous emission. The inverse dependence of the relative intensity fluctuations 
on the intensity is characteristic of the optical laser. It is expected to be true 
for the atom flux from atom lasers. In contrast, for thermal beams the relative 
intensity fluctuations are constant. 

The phase of an optical laser far above threshold is also stabilised. 
This is measured by the electric field amplitude temporal coherence func- 
tion {E* {ti)E{t2)) ■ Its Eourier transform is the spectral density, which is 
approximately LorentziancJ. Its width is the_Laser bandwidth, which is found 
to be inversely proportional to the intensity^. This is another characteristic 
optical laser property which we would expect to be true for atom lasers. 

We now describe the coherence functions which are used to quantify the 
coherence of fields. The amplitude coherence function is: 

(A*(xi)A(x2)), (17) 

where the angle brackets denote an average over noise, and a bold denotes 
the spacetime event x/j = {xk,tk)- This coherence function is unambiguously 
referred to as a two-time amplitude-correlation functimi. In the literature it 
may be referred to as either a firstO or second ordeiO amplitude coherence 
function. We use the former convention. If this function factorises 

(A*(xi)A(x2)) = (A*(xi)) (A(x2)) (18) 

the amplitude is said to be first order amplitude coherent. This definition 
captures the concept that a coherent field has no more noise in amplitude 
products than that due to the noise in the amplitude itself. 

The definition of coherence as the degree to which coherence functions 
factorise is retained for quantum fields. Coherence furictions (also called cor- 
relation functions) of order n were defined by Glaubertj in terms of the field 
operators by 

G'("Hxi . . . , X„, X„+l . . . X2„) = (*^(xi) . . . *^(x„)^'(x„+l) . . . *(X2„)), 

(19) 

where the angle brackets denote the quantum mechanical average. If G(") 
factorises 

G(")(xi . . . ,X2„) = (vl/t(xi)) . . . (vl/t(x„))(vl/(x„+i)) . . . (vl/(x„)) (20) 

the field is said to be nth order coherent. This particular form of the coherence 
functions arose from photodetection theory. The pairing of creation and an- 
nihilation operators is a result of practical photodetectors measuring photon 
number. The normal ordering, in which all creation operators are to the left of 
all annihilation operators, arises because detection relies on absorption rather 
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than stimulated emission, which would be affected by spontaneous emission. 
An all orders coherent quantum state is the annihilation operator eigenstate, 
the so called "coherent state". This is the quantum state that most closely 
approximates (superpositions of) noiseless classical waves. 

One difficulty with this definition of coherence is that it relies on the 
spontaneous symmetry breaking assumption that the field expectation (^'(x)) 
value is not zero. Normalised coherence functions overcome this problem: 

(n), ^ G(")(xi...,X2„) 

^(xi . . . ,X2„) = — ^ — -jj^. (21) 



(G(1)(X1,X1)...G(1)(X2„,X2„)) 



For a maximally coherent field the numerator factorises and g^"> — 1. 

In a simple spatial interference experiment two sources at xi and X2 
propagate for the same time to a detection point. If the sources are co- 
herent, interference results. The relevant coherence function is first-order: 
g^^^xi, X2), because the intensity at the detection point contains the cross 
term {^''' (xi)^ {X2)) in the field amplitudescJ. (Note that when either the 
spatial or temporal arguments are suppressed, as are the temporal arguments 
here, they are assumed to be equal.) If this cross term is zero the field fluc- 
tuations wash out the interference. The visibility of the interference fringes is 
maximal when g^^''{xi,X2) = 1. From the definition Eq.(^l|), this is true for 
any single mode field, as for (7'^^ the only relevant expectation value is {a^a). 
In a BEG a large number of atoms are cooled into a single mode. Therefore a 
BEG with output coupling that preserves coherence is a prime candidate for 
an atom laser. 

The first order spatial coherence of two freely expanding Bose-Einstein 
condensates has been observed experimentally by the detection of interference 
fringestZI. In this experiment two condensates separated by 40/im were created 
in a double- well potential. The potential was switched off and the condensates 
left to expand over a period of 40ms after which time they overlapped. High 
contrast matter-wave interference fringes were observed indicating that the 
BEGs were first order spatially coherent. 

Ideally one would want a continuous output beam which is stationary, i.e. 
for which the first order temporal coherence, g''^\t -f r, i) = ^''^■'(t) does not 
depend on t. The first order temporal coherence function for an atom laser 
beam, at some particular point, is defined as 

:^^^'r)- , + (22) 

This is a useful measure of phase fluctuations, which we require to be small 
for a laser like source. For a laserlike, first order temporally coherent source. 
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g^-^'^ir) « 1 over a time called the coherence time. This is in contrast to a. 
noisy source for which g^^^r), decays quickly to zero. According to Wisemanlll 
an atom beam might be regarded as laser-like if the coherence time is much 
greater than the inverse atomic flux. 

The second order coherence function relates two separate field quantum 
(atoms or photons) detection events. Unlike the first order coherence function, 
it is sensitive to the quantum mechanical state of the field. For a single mode 
thermal state g^"^^ = 2 and for a BEC or single mode coherent state g^"^^ = 1. 

The second order spatial coherence of a matter atave source is related to 
the atom-atom interaction energy of the condensatecj. This is because for a 
short-range potential, the interaction energy is proportional to the probability 
that two atoms ate nearby, which is in turn proportional to g^^^ . The exper- 
imental evidenceEJ is consistent with g^^' = 1.0 ± 0.2 for BECs, confirming 
that condensates suppress local density fluctuations. 

As well as first order temporal coherence, an atom laser beam is expected 
to have higher order temporal coherence. The second order temporal coher- 
ence function (for a stationary system) is 

where / = ^^5', and the colons denote normal ordering, i.e. all creation oper- 
ators to the left of all annihilation operators, g*^^^ is approximately unity for 
a coherent source. A filtered thermal beam, which is first order coherent, has 
^(^^(t) = 2 for short times, approaching 1 for long times. In terms of counting 
field quanta g^^^r) is determined by the distribution of arrival times. A field 
with ^'■^''(t) — 1 has arrival times distributed according to Poissonian statis- 
tics, while g^^^r) > 1 corresponds to superpoissonian statistics, or bunching 
of arrival times. Only purely quantum mechanical correlations, such as those 
associated with squeezed states, allow ^'^'(t) to be less than one. 

Experimental evidence of third order coherence in BECs has been given by 
Burt et a/.Ell. They compared the three-body recombination rate constant in 
condensed and non-condensed Bose gases. They found that the ratio of g'^^^ for 
a therpial cloud to that for a BEC was 7.4±2, in agreement with the predicted 
value E3 of 3! = 6. For a BEC .g^^) 1, and for the non-condensed fraction 
g^^") = 6. It is notable that g^^'^ = 1 is the same as for classical distinguishable 
particles. This is because BEC particles occupy the same field mode and are 
therefore already indistinguishable in principle. Hence state symmetrization, 
which is used to enforce indistinguishability, thereby generating the quantum 
mechanical factor of 3!, is neither necessary nor possible. 

For light only the Glauber type coherence functions that we have discussed 
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are relevant, because all light detectors use photon absorption. However atoms 
may be detected in a variety of ways, and different coherence functions may be 
measuredcJcJ. For example the spectrum of off resonant light scattered from 
a BEC is related to the density correlation function, rather than the Glauber 
normally ordered coherence functionES. Normally ordered coherence functions 
are relevant for destructive detection, such as by ionisation detectors. We have 
focussed on them because they allow the most straightforward comparison 
with the optical laser. 

5 Atom Laser Models 

In this section we discuss the types of theoretical models which have been 
used to describe the atom laser. The starting point is usually an analogy with 
the optical laser. However despite the similarities, the differences between 
atom and optical lasers are also of considerable interest. We consider the 
advantages and disadvantages of three major types of atom laser models. In 
order of increasing complexity: rate equation, mean field, and fully quantum 
mechanical. 

Rate equation models only consider state populations, ignoring the dy- 
namics of coherences between quantum states. Hence they are unable to 
properly describe the quantum mechanical coherence properties of the atom 
laser. Nevertheless they usefully describe important physics such as the role 
of Bose-enhancement in producing the laser threshold. 

The mean- field approach uses Gross-Pitaevskii type equations, Eq. (|l5|) . 
Since the atom laser is then described by a macroscopic wavefunction, quan- 
tum coherences may be included in the models. Additional terms may be 
introduced to describe the pumping and output coupling. A limitation of 
these models is that they assume that the field is in a coherent state, as dis- 
cussed at the end of section ^. This means that the quantum coherence of the 
output cannot be calculated. However the nonlinear dynamics of the BEC and 
output beam arc fully modelled, as is the spatial propagation of the output 
beam. 

The most common fuJly quantum mechanical approach uses the quantum 
optical master equationE3. This is an equation for a reduced density operator, 
obtained by averaging (strictly tracing) over part of the system, usually the 
output and/or pump modes. The Born-Markov approximation is commonly 
applied during the trace over the output modes to ensure that a_differential 
equation, rather than an integro-differential equation, is obtainecS. Because 
this approximation is not necessarily valid we shall also discuss more general 
quantum operator models. 
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Figure 2. Schematic level diagram for the atom laser model of Moy et a/J2j (left). Plot of 
the steady state number of atoms in the lasing mode, as a function of the dimensionless 
pumping rate, R (right). Threshold occurs at iJ = 1. 



In the rest of this section we give brief descriptions of these three ap- 
proaches and the kinds of physics they treat. Pumping and output couphng 
are considered in detail in sections 6 and 7. 

5.1 Rate equations 

Thi! first rate equationJpased atom laser models were presented by Olshanii et 
al.U and Spreeuw et alu. We discuss a later model due to Moy et aZ.E3 which 
is based on implementation in a hollow optical fibre. 

The model consists of atoms with four energy levels, as indicated in Fig. ^. 
Level |1) is the input pump level, level |2) is the lasing level and level |4) is 
the output level. Level |3) mediates the output coupling Raman transition 
and is not populated. Level |1) atoms irreversibly emit a photon to populate 
the lasing level |2) at the Bose-enhanced rate ri2(iV2i + 1), where is the 
number of atoms in level |2) and the ground state of the laser trap. A Raman 
transition couples atoms from level |2) into level |4) with rate constant r24- 
The resulting rate equations for the numbers of atoms Nk in level \k) are 



dNi 



Y,gj{N2, + l)ri2iVi - (1 -Y^g^y^^N, 



dt 




dN2j 



9jri2Ni{N2j + l)~r24N2j 



+ G,{N2j + l)N4r24, 



dt 
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^ = E ^2^-^24 - E Go{N2, + l)iV4r24 " N^^. (24) 

The iV2j are the numbers of atoms in level |2) in the jth mode of the lasing 
level trap, gj is the overlap between the pump mode and the jth mode of the 
lasing trap. Gj = N^j/Ni is the proportion of atoms with electronic level |4) 
that have the same spatial mode as the jth mode of the lasing trap, to is the 
time scale on which atoms are irreversibly lost from the system due to the 
Raman momentunp_kick. For a more complete description of these equations 
refer to Moy et aLE£. 

The solutions of these rate equations show behaviour familiar from optical 
lasers. The steady state population of the laser mode N21 as a function of 
a suitably scaled pump rate R = 7'igi/?'24 shows a rapid increase at the 
threshold value of i? = 1, as shown in Fig. |^. As a function of time the 
above threshold values of the lasing cavity mode populations N2j show initial 
growth, and then all but the lasing mode decay. This is the familiar mode 
competition in which the "winner takes all" dueJ;o runaway Bose stimulation. 

The rate equation model of Olshanii et ala is similar in form. However 
they included photon reabsorption and found that lasing is inhibited by small 
amounts of reabsorption. 



5.2 Mean- field models 

In this approach, the atom laser is modelled using the Gross-Pitaevskii (GP) 
type equations discussed in section |[ The "mean-field" is the GP wavefunc- 
tion which determines the effective potential felt by an atom due to its 2-body 
interactions with all of the other atoms in the condensate. Indeed one of the 
advantages of GP models is this straightforward treatment of atom-atom in- 
teractions. _ 

Kneer et aLc3 describe a GP model of an atom laser including pump 
and loss terms. The modified GP equation for the laser mode macroscopic 
wavefunction is 

= + VEir)4' + C/olV^lV + i^gain0 + ^lossV-- (25) 

The first three terms are standard. The gain and loss are represented by 
phenomenological terms with rate constants F and 7c respectively, 

Hgain^ = yFA^^V, ■H'lossV' = ^y^cV'- (26) 
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In addition to the GP equation (25), a rate equation is used to describe the 
number of uncondensed atoms, N^^ 

^ = - luNu - TNcNu. (27) 

Ru describes a source which pumps the uncondensed state at rate Ru- The 
term juNu describes atoms lost from the system, but not trapped in the 
condensed state. The final term describes the transfer of atoms into the con- 
densed state. The condensate loss is spatially homogeneous, and undamped 
collective excitations are found. To avoid this they introduce a spatially de- 
pendent decay rate, so that 7c is a function of position, 7c (a;). This produces 
a single lasing nmd^ as the spatially dependent loss leads to mode selection. 

Robins et al£3 present a GP model which combines the pumping model of 
Kneer et al., with an explicit treatment of the output coupling and of the out- 
put beam. Including the output beam in the model allows one to investigate 
its spectrum and spatial distribution. Furthermore, non-Markovian effects in 
the output, such as Rabi oscillations between the condensate and beam, can 
be explored. The equations which describe the laser mode ip°'{x,t), and the 
output beam ip^{x,t) are 

-ilr{\r\^ + IV-'H^' - ?ir^e-^"°^V°. (28) 

The number of uncondensed atoms obeys Eq.(p7|). The Raman output 
coupling is described by the terms containing the momentum kick generators 
^±ikox ^ Three body recombination has been included as a condensate loss 
mechanism, by the terms proportional to 7^. This is found to have a stabilising 
effect on the non-linear dynamics of the atom laser. Two-body losses have 
not been included. 

Typical output from a one dimensional numerical solution of this model is 
shown in Fig. ^ The oscillations represent nonlinear dynamical noise, which 
will determine the laser's effective linewidth, or first order coherence. This is 
an example of physics which cannot be treated by the rate equation approach. 

How mean field models may be extended to incorporate finite temperature 
effects, such as the population of excited levels, is described in subsequent 
sections, and by Burnett and Griffin in this volume. 
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Figure 3. Atom laser beam density as a function of position and time, showing a quasi- 
steady statei — lAtoms are absorbed at the spatial boundaries. Parameters are given by 
Robins et aijl3. 



5.3 Fully quantum mechanical models 

A number of quantum optical masteLequaJiDn based models of the atom laser 
have been discussed in the literatureOO'Ht3. The model of Holland et alu is 
based on evaporative cooling and incorporates: the creation of an atom in the 
trap from the pump field, the loss of an atom from a high lying trap state, 
collisions which scatter atoms into or out of the ground state, and the loss of 
atoms from the ground state through output coupling. 

A three level model is sufficient to describe these features. Atoms are 
injected with rate ki in level one and are lost with rates kq and K2 from level 
zero and level two respectively. Only the energy conserving 2-body interaction 
terms are considered. Level two is adiabatically eliminated by assuming large 
damping due to the evaporative cooling. This leads to the master equation 
for the density operator p, 

% = ^\^^^P\ + ^o^^\a^]p + ^^Mo^i]p + ^Malafip (29) 

where 

= h{alao - 2no)(Vboooa(')ao + Voioiajai) -I- hViiiiafa\, (30) 

and no is the mean population of level zero. The loss terms are assumed to 
be of the standard (Born-Markovian) Lindblad form, defined by 

V[c]p = 2cpc'' - {c^cp + pc^c) (31) 
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The effective redistribution rate fir is proportional to |Vo2iiP/k2. 

Equations for the populations, {nk\p(t)\nk), where \nk) describes a state 
with n atoms in level zero and fc = or fc = 1 atoms in level one, are 
obtained. Level one is rapidly depleted by two atoms scattering into levels 
zero and two. The steady state population of the level zero lasing level is 
found to be Poissonian with mean atom number 2ki/(3ko)- Thus when loss 
out of the lasing mode is sufficiently small kq << ki, a large population can 
build up in the lasing mode. 

In order to produce a state with a well defined phase they consider pa- 
rameter regimes where J7r is much larger than the frequency shifts in Vi . This 
inhibits the dispersive nature of the latter from destroying phase coherence. 
They find that, as for an optical laser, line narrowing occurs as the population 
of the ground state increases. However, in a refinement of the model due to 
Wiseman et aZHj the linewidth of the laser output is broader than the bare 
linewidth of the laser mode. 

The crucial refinement is in the description of the pumping of level 1 . The 
approximation that Holland et al. make allows atoms to leak into, but not 
out of, the trap. To allow the latter Wiseman et al. replace the term KiI?[a|]/9 
with 

Ki {N + l)V[a\]p + KiNV[ai]p. (32) 

The N is the mean population of the pump reservoir modes. In fact, there is 
no limit in which the extra terms can be ignored. When kq << ki, there are 
two distinct parameter regimes depending on the relative size of a parameter 
proportional to K2|Vb2iiP, which they call F, and kq. These they call the 
strong and weak collision regimes: T » kq and T << kq respectively. In 
the weak collision regime, T << kq there is a large population in the laser 
mode if N » 1, i.e. with strong pumping. The power spectrum reveals a 
linewidth larger than the bare linewidth of the lasing mode. Nevertheless, the 
phase diffusion rate may still be slow in the sense of being much less than the 
output flux of atoms frorn-the lasing mode. 

Wiseman and CoUetto proposed an atom laser model using optical dark 
state cooling. The laser is found to have a linewidth narrower than the bare 
linewidth. This is an improvement on models based on evaporative cooling 
discussed above. The output is also found to be second order coherent. 

The master equation is one means of describing the full quantum statis- 
tics of an atom laser. It has the great advantage of drawiiLg^on the wealth 
of techniques that have been developed in quantum opticsEj. However, a 
fully quantum theory of the atom laser caii J:ie_developed without using the 
standard master equation approximationsCj^Ell. In particular, assuming a 
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Born-Markov form for the output coupling ignores the possibihty of the out- 
put beam acting back on the laser mode. This may happen, for example, due 
to Rabi ciZEling between the output and laser modes. 

Hopeo developed a general theory of the input and output of atoms from 
an atomic trap. Moy et al£3 used this theory to find an analytic expression for 
the spectrum of a beam of atoms output from a single mode atomic cavity. For 
very strong coupling it is just the spectrum of the trapped lasing mode, and 
is therefore broad. As the strength of the coupling is reduced, the linewidth is 
reduced. For small coupling strengths the lineshape approaches a Lorentzian, 
centred about the energy of thft_cavity with a linewidth proportional to the 
coupling strength. Hope et alEU include a pumping mechanism and they 
find that the output spectrum is substantially wider than that obtained using 
the Born-Markov approximation. However an important limitation of these 
theories is that they do not include 2-body interactions. 



6 Atomic Gain 



We have seen in sections |^ and ^ that a mechanism for atomic gain, to provide 
coherent amplification of the atomic field in the lasing mode, is an essential 
element of an atom laser. A number of gain mechanisms have been proposed, 



(e.g. see section 5.1) but to date only one, evaporative cooling, has been 
demonstrated experimentally. The process of evaporative cooling is based on 
the removal of high energy atoms from the gas, so that subsequent collisional 
rethermalisation of the remaining atoms results in a distribution of lower tem- 
perature (i.e. lower average energy). Under appropriate conditions, bosonic 
stimulation of transitions into the ground state of the confining potential leads 
to the formation of a Bose condensate: in an atom laser this would constitute 
the lasing mode. 

The kinetics of condensate growth has long been a subject of theoreti- 
cal study^^~^^, but the first quantitative calculations were carried out only 
recently by Gardiner. Together with coworkers he has given a comprehen- 
sive treatment of condensate formation by bosonic stimulation, in a series of 
papers^*^"^^. In this section we review the major features and results of this 
theory. 



6.1 Quantum Kinetic theory 

The starting point of the treatment is the third in a sejdfis of papers on Quan- 
tum Kinetic theory developed by Gardiner and ZoUeiO, which we shall refer 
to as QKIII. We note that an alternative formulation of a Quantum Kinetic 
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Figure 4. Schematic representation of modification of trap energy levels by condensate 
meanfield. Levels above En are unaffected. 



theory has been given by GrifRn and coworkers (see this volume, and refer- 
ences therein). The gas of Bose atoms in a harmonic trap of frequency lu is 
described by a second-quantized field, as in Eqs.(]l| - 14|). Two-body collisions 
are included but three-body collisions are neglected. Guided by the fact that 
the energy spectrum above a certain value En is unaffected by the presence of 
the condensate, we can divide the field into the condensate band (all energy 
levels below Eji) and the non-condensate band (all the energy levels above 



The noncondensate band (which contains the majority of the atoms) is 
taken to be fully thermalised, in the sense that for any point r, k in phase 
space the single particle distribution function F{r, k) (in principle a Wigner 
function) is localised, and has a well defined temperature T and chemical 
potential /i. Thus the noncondensate band can be treated as a time inde- 
pendent bath, and a master equation is derived for the density matrix of 
the condensate band, (Eqs. (50a)-(50f) QKIII). It is worth emphasizing that 
in general the condensate band is not in thermal equilibrium, and the mas- 
ter equation provides a description of its dynamical evolution. The level 
structure of the condensate band is well described by the Bogoliubov spec- 
trum, and GardinerH has shown that appropriate states for this band are 
l-^j {i^m}) which have a definite number N of total particles in the band, and 
a set of numbers {um} representing the occupation numbers of the conden- 
sate (m — 0) and quasiparticle (excited Bogoliubov) levels. In this treatment, 
the condensate wavefunction and energy eigenvalue (the condensate chemical 
potential /ic(-/V) ) are given by the time-independent Gross-Pitaevskii equa- 
tion with hq atoms. The quasiparticles may be particle-like excitations (at 
higher energies) or phonon-like collective excitations (at low energies). The 
phonon excitations are not eigenstates of the number operator and thus the 
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Eji), as represented in Fig. ^. 



total population of the quasiparticles need not be conserved, in contrast to 
iV which is strictly conserved. 

By taking a representation of the master equation in terms of the states 
\N, {rim}), and then taking the diagonal elements, a stochastic master equa- 
tion for the occupation probabilities p{N, {um}) is obtained. The transitions 
appearing in this equation can be divided into those which change the num- 
ber of atoms in the condensate band, and those which do not. The former, 
which we call growth processes, arise predominantly from collisions of two non- 
condensate atoms where one of the outgoing atoms goes into the condensate 
band. The reverse collision is also included in the growth processes. The 
transitions which do not change N are called scattering processes, and arise 
predominantly from collisions between a condensate band atom and a non- 
condensate band atom which leave one atom in each band. Other collision 
processes, in which the pair of colliding atoms either begins or ends in the 
condensate band, will in principle also contribute to growth or scattering, but 
can be neglected since the condensate band contains such a small fraction of 
all the atoms. We note that collisions between pairs of atoms that begin and 
remain in the noncondensate band play the role of thermalising the bath and 
do not appear in this master equation. 



6.2 'Particle-like' rate equations 

Using standard methods, a deterministic rate equation for the mean popu- 
lation {um) of the mth quasiparticle level (which we will subsequently write 
Um) can be obtained from the master equation and we write it as 

— '^mlgrowth + '^mlscatt ("^"^) 

where the forms for the growth transition rates nm| growth and the scattering 
transition rates nmlscatt will be discussed below. A similar deterministic rate 
equation can be written for the mean population of the condensate band (N) 
(which we will subsequently write N). 

The growth transition rates growth can be obtained in explicit form 
from the population master equation given in Eq.(189) QKIII, and comprise 
all the processes where TV changes by 1, and one of the quasiparticle popula- 
tions, Uj say, changes by 1 or 0, (with all other rife unchanged). Most of these 
processes (e.g. N ^ N + 1 and Uj rij + 1, or N ^ N + 1 and {um} un- 
changed) admit a particle-like interpretation. However, those processes where 
N and Uj change in opposite directions have to be interpreted as phonon-like 
processes, and for example could represent a process where the condensate 
is initially oscillating, and a thermal atom is added to the condensate band 
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and causes the condensate oscillation to be reduced. All such phonon-like 
processes are neglected in this treatment, since only a small fraction of quasi- 
particle wavefunctions have phonon-like character, and we expect they will 
play an unimportant part in the condensate growth. The growth terms for 
the quasiparticle levels thus simplify to the form 

Igrowth = 2VK++ (N) (n™ + 1) - 2W--iN)nrn (34) 
and the equation for the mean number in the condensate band becomes 
dN 

— = 2W+{N){N + 1) -2W-{N)N (35) 

where W^~^{N) and W'^{N) are forward rates and W~~{N) and W~{N) are 
backward rates. It is shown in QKIII that when the condensate wavefunction 
fAr(r) is sharply peaked compared to the spatial dependence of F(r,k), then 

W+iN) - ^2^^^^3 J '^'K:! / '^'Ks J rf^Ka J d^k6{Ki + K2 - K3 - k) 
S{Au;u3 " fic{N)/h)F{0, K,)F{0, K,) (1 + F(0, K3)) ||^(k)|2 (36) 



The integrand represents a collision in which a pair of noncondensate atoms 
of momenta Ki and K2 produces one noncondensate atom (K3) and one 
condensate band atom of momentum k. Energy conservation is expressed in 
the second 6 function, in which Aa;i23 = "-^Ki + ^ "^k.,, where Tiuj-k- = 
?i^Kf/2m +Vt(0) is the particle energy at the centre of the trap potential Vr- 
^iv(k) is the momentum-space ground-state wavefunction. The expression for 
Wj^"*" (N) is similar to Eq. , and the forward and backward rates are related 
by 

W+{N) = e[''-''c.Wl/'='^W^-(7V) (37) 
W++iN)^ei'^~'-^/''^W,--iN) (38) 

where e™ is the energy of the m*'' quasiparticle level. Equations (^3|)-(p8|) 
are derived within the approximation that the number of particles in the con- 
densate, no, is large enough that we can write uq « N, (which is required 
for the Bogoliubov spectrum to be accurate). However, in the early stages 
of condensate growth when uq is small, the unperturbed spectrum should be 
used, and will be accurate provided none of the condensate band levels are 
too highly populated. This means that by making the replacement N hq in 
Eqs.(p3|)-(^8|) those equations will provide a good description of the dynamical 
evolution of the levels in the condensate band throughout the formation of 



Ballagh: submitted to World Scientific on February 1, 2008 20 



the condensate. Values for the rates can be obtained in varying levels of 
approximation, and_the simplest expression, used in the first quantitative de- 
scription of growthtS, is obtained by assuming a classical Maxwell-Boltzmann 
distribution for F(r, k), and allowing the integrals in K.1,2,3 in Eq.(|6|) to ex- 
tend over all energies, rather than just the noncondensate band. This leads 
to the explicit expression 



4m(afcr)2^2M//=T 



(39) 



in which we note the prefactor to the exponential is independent of uq and 
is essentially the elastic collision rate pav with quantities_eiialuated at the 
critical point for condensation. More accurate evaluationsoEj of W'^ which 
take into account the Bose-Einstein nature of the non-condensate distribution 
and exclude the condensate region in the K integrals, give a value about three 
times larger than that of Eq.(^9|). iy++(no), which has a form similar to 
Eq. ( ^6| ) , is more difficult to evaluate but is of the same order of magnitude as 
W^{no), and in the first instance we set lV,++(no) « M^+(no). The sensitivity 
of the results to this approximation will be discussed in section 3.2. Finally, 
in order to reduce the number of equations to a tractable set, an ergodic 
assumption is made which allows us to group together levels of similar energy 
into small sub-bands, with only the ground state being described by a single 
level. The sub-bands have (central) energy Cm and degeneracy gm, and the 
equations for the growth processes of these bands are 



"mlgrowth = 2iy"++(rio) I 

"olgrowth = 211/+(no) I 



1 _ g[e,„-M]/fcT' 



1 _ g[^^c{no)-^i]/kT 



no + l} 



(40) 
(41) 



These equations exhibit the most important features of condensate growth, or 



atomic gain, as we shall discuss shortly in section 5.3. The remaining terms 
required for Eq.(^3|), the scattering terms jimlscatt: describe processes in which 
a condensate band atom and a noncondensate band atom collide and leave 
one atom in each band. Gardiner and LeeEj have evaluated the scattering 
terms within well defined approximations (for the case of an isotropic trap), 
to give 



Iscatt 
k k<m 



ttTi 
nk{g 



)e-''"™'=/'=-^-n,„(5,+nfc) 
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— \nki9m + rim) - n^igk + nk)e ^^-^/^'bT 



kyrn 



9k 



(42) 



where T{T) — e '^"•/'^'^ , fiuJnq = e„ — Cg, and the notation k > 



means ek > Sr, 



6.3 Results 

Equations - ^) provide an exphcit form for the set of deterministic rate 
equations for populations in the condensate band, with the size of the set 
determined by the amount of level grouping. Much of the behaviour that 
occurs in the full set of equations is captured by a simple imnlementation 
in which the condensate band is represented by only one leveEa, i.e. the 
condensate. The equations then reduce to a single equation (p|), which has 
the laser-like features of stimulated transition rates (the terms proportional 
to no) and a spontaneous transition rate (the term +1 in the curly braces) 
which initiates growth. In addition, this so-called simple growth equation has 
thermodynamic features, in that the forward and backward transition rates 
become equal (i.e. equilibrium is reached) when the chemical potential of 
the condensate, /ic(no), becomes equal to /i of the thermal bath (to order 
1/N). A crucial aspect of the growth process is that as the condensate grows, 
the meanfield causes /ic(rio) to increase, beginning from /ic(0) = 3hLu/2 (the 
trap zero-point energy), and rising as (no)^^^ (in the Thomas- Fermi approxi- 
mation) at large no. It is clear that a condensate can only grow if the thermal 
bath has /i > 0, and this is achieved in the process of evaporative cooling by 
removing the highest energy atoms. The higher energy levels quickly come 
to their equilibrium distributions, but the lower energy levels remain far from 
equilibrium, and the ensuing tendency of collisions of lower energy bath atoms 
to feed the condensate band can be characterised by a bath chemical potential 
^l>0. 

A typical result for the simple growth equation is shown as the dotted 
line in Fig.^ where we have chosen a case corresponding to the conditions 
of the MIT growth experimentEZI in which sodium atoms of scattering length 
a = 2.75 nm are contained in a trap of geometrical mean frequency a; = 27rx 50 
Hz. Although we do not expect quantitative accuracy from this simple model, 
the characteristic S shape of the growth curve exhibits the key aspects of the 
growth process. In the initial stages the growth is slow and arises from spon- 
taneous transitions; then as no increases, stimulated transitions dominate and 
the growth is steep. Finally, as the condensate chemical potential approaches 
fj,, there is a slow approach to equilibrium. 
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Figure 5. Typical condensate growth curves for different implementations of the model; 
simple growth equation (dotted line); full model (solid line); with no scattering (dashed 
line). Parameters are appropriate to the regime of the MIT experiment: = 43.37ia;, 
T = 900 nK, no{t = 0) = 100. 



Of course for a quantitatively accurate treatment, more condensate band 
levels must_be included, along with the more careful calculation for W'^ (see 
Eq.(lO) reiEj). Typically 20-50 sub-bands are used, each of width « Ttu. The 
energy e,n of these sub-bands increase as no grows, with the lower levels most 
affected, and the levels at the top of the condensate band not affected at 
all. Although the increase of the condensate energy level is crucial to the 
growth process and must be treated as accurately as possible, it is found that 
there is little sensitivity to the scheme used to alter the quasiparticle energy 
levels. Typically, the levels between fj,c{no) and 2fic{n-o) are compressed 
from the bottom in a linear fashion, and higher levels are unaltered. The 
solid line in Fig.|| shows the effect of using a set of sub-bands in the growth 
process. The physical parameters are identical to those used in the simple 
growth model (the dotted curve) but there has been a speed up of about 
an order of magnitude. The effects of the scattering terms can be separated 
from the growth terms, as we show with the dashed line in Figj5| which is a 
calculation using the same number of sub-bands as for the solid curve, but 
with all scattering terms put to zero. It is clear that the main effect of the 
scattering processes is to cause the initiation of condensate growth to occur 
much more sharply. This is due to the fact that the growth processes can 
initially populate a wide number of sub-bands, and then scattering allows 
these populations to be redistributed directly into the condensate. The steep 
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part of the growth curve is essentially unchanged by scattering, since it is 
dominated by the stimulated terms in the growth processes. |— . 

A number of sensitivity studies were carried out by Gardiner and LeeEj, 
who showed that the scattering transition rates could be varied by two orders 
of magnitude, the relative size of W'^~^ to could be varied by a factor 
of 2, and a wide variety of initial conditions for the rate equations could be 
used, with only moderate effects on the growth curves. 

6.4 Full dynamical solution 

The master equation approach described above, which assumes a nonde- 
pletable bath of fixed thermodynamic properties, can be expected to be in- 
accurate when high fraction condensates form. Under those circunistances 
the full dynamics of the system must be considered, and Davis et al. have 
treated this problem using a modified Quantum Boltzmann formulation. In 
their treatment the single particle distribution function /„ for the n*'* energy 
sub-band evolves according to 

9n — "^^^ / ^ "e,,+e„,ep+egS'min(p,(j,r,n) 

P-Xi-X 

X {/p/, (1 + !r) (1 + /„) - (1 + /p) (1 + /,) /r/4 . (43) 

This equation appears similar to the ergodic form of the well-known Quantum 
Boltzmann equation (e.g. see ref Ej) but the crucial modification from Quan- 
tum Kinetic theory is that the energy levels Cm and the degeneracies alter 
as the condensate population changes. A division of the levels into condensate 
and noncondensate bands is now made only for computational purposes, and 
the evolution of all sub-band levels is explicitly followed. In addition to the 
collisional processes treated in the master equation method, Eq.(^3|) includes 
thermal scattering processes where the pair of colliding atoms begin and re- 
main in the noncondensate band, as well as collisions which either begin or 
end with a pair of atoms in the condensate band. As before, all levels are 
treated as particle-like, which allows an analytic expression to be obtained 
for the density of states and its dependence on the condensate occupation. 
The numerical simulation of Eq.(^), which we shall refer to as the ]ull dy- 
namical solution, requires vastly more computation than the rate equation 



model described in section 6.2, and considerable care must be taken with en- 
suring energy conservation in the collisional calculatians and in rebinning the 
sub-bands at each time step. Details are given in ref Ej. 

For small condensate fractions (of less than about 10%) the full dynam- 
ical solution agrees well with the rate equation solution. However for larger 
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Figure 6. (a) Comparison of condsDsate growth models for a final condensate fraction of 
24%. Solid line, simulation of Eq. (43) [with /i = before evaporative cut]. Dotted line, rate 
equation model. Final condensate parameters ng = 7.5 X 10^, T = 590nK. (b) Chemical 
potentials ^cff' (solid line) and /ic("o) (dashed line). 



fractions, as illustrated in Fig. ^(a) we find that the bath dynamics have an 
appreciable effect. For this case, where parameters are chosen to be appro- 
priate for the MIT experiment and the final condensate fraction is 24%, the 
full dynamical solution (solid line) grows faster than the rate equation solu- 
tion (dotted line). This can be most easily understood by introducing the 
concept of an effective chemical potential /icff for the thermal cloud (obtained 
by fitting a Bose-Einstein distribution to the distribution function above Efi) 
and using this in the simple growth equation Eq.(pl|). Figure ^b) shows the 
evolution of both fies and /ic("-o) during the simulation, and we can see that 
during the stimulation-dominated period of growth, /icff substantially exceeds 
the final equilibrium value, and thus from Eq.(|4]]), the backward rate from 
condensate to thermal cloud is greatly inhibited. Ultimately this is due to the 
very severe cut that is required on the initial equilibrium distribution in order 
to allow evaporative cooling_to proceed to such a high condensate fraction. 
We note that Zaremba et aZ£il have treated the full dynamical problem using 
a somewhat different formalism, but obtain essentially the same numerical 
results. 

The theory reviewed in this section |CaxiL.h£j compared to the experimental 
results of the MIT groupLll It is found EjCJ'Ej that there is very good agree- 
ment in some parameter regimes, but significant discrepancies in others, which 
may be attributable to the difficulty of fully characterising the experimental 
parameters. A definitive comparison is handicapped by the unavailability of 
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a wide range of primary experimental data, and further experiments will be 
required to comprehensively test the theoretical predictions. 



7 Output Coupling 

An output coupler for an atom laser is a device which extracts atoms from the 
laser mode while preserving their coherence. All current coupler mechanisms 
are based on applying an electromagnetic field to-change the internal state 
of the atom from a trapped to an unttapped statm. The first demonstration 
of an atom laser, by the MIT groupE3, used an RF pulse to transfer atoms 
into untrapped zecman substatcs, where they then fell freely under gravity. 
A refinement of this technique, using a precisely controlled trap and_cw-RF 
fields has been used to produce a quasi-cw coherent output beamEJ. The 
use of Raman coupled output fields offers some advantages, principally in 
providing a momentum kick to the outgoing atom beamE^I, but also in that 
the optical fields can be focussed into a small volume, achieviag a measure of 
spatial selectivity for the coupling. Recently the NIST groupEZI demonstrated 
Raman output coupling in an experiment where the output beam was ejected 
in the horizontal direction. i— in 

A number of meaiL field theoreticaLtxeatments have been given of RFc2lo 

et 



and Raman couplersl__l, and elsewheii 
eluded. Recently the Oxford groupEZI 



the effects of gravity have been in- 
3, have given a detailed treatment of 
output coupling from a finite temperature condensate. In this section we re- 
view the theory of output couplers, beginning with a mean field treatment, 
which establishes a framework for both RF and Raman output couplers. We 
then briefly summarize the theory of finite temperature couplers. 

7. 1 Mean field treatment of output coupling 

The basis of the mean field treatment of an output coupler is a pair of coupled 
GP equations for the mean fields ipi and 'ip2 of atoms in internal states |I) 



and |2), as discussed previously in section 5.2. For convenience, we present a 
dimensionless form of the equations 

^=iV2V'i-i!/iVi+i^^Y^^2-iC[|V^iP+«;|V2|']V'i-*G#i, (44) 

^=^S/^^2-iV2iJ2+i^^^^l-^C[\tP2\^+w\i:l\^]iJ2-^Gy^J2 . (45) 

We shall assume that state |1) is the trapped state, with a harmonic potential 
Vi while state |2) is untrapped (V2 — 0) or anti-trapped (V2 = —Vi), so 
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that '02 represents the output field. In Eqs.(}44D and ( [45[ ) the units of time 
(to) and distance (io) are respectively and {h/2mujy/^ where to is the 
trap frequency, and the fields are normalised such that /[|V'iP + l'02p] d'^'T = 
1. The nonlinearity constant C is 4:TThaN/{mujXQ) (in 3D), and the relative 
scattering length w of the inter- to intra-species collisions is typically close to 1. 
Gravity acts in the —y direction, and the scaled acceleration is G = mgxo/hw. 
The coupling term A can describe either a single photon coupling or a two 
photon (Raman) coupling, and can be written in terms of a slowly varying 
envelope r2(r, t) 

A(r, t) = n{r, t) exp i(ke,„ • r- Agmt) (46) 

where kem and Aem measure the net momentum and net energy transfer 
from the EM field to the output atom (in units Ti/xq and Tiuj respectively). 
For the case of a single photon transition, Aem = loi — uj2i, where uji is the 
photon frequency, W21 is the energy difference E2 — Ei between states |2) and 
|1) (in the absence of the traps), and ri(r,t) is the Rabi frequency. For a 
Raman transition connecting |1) to |2) (through a far detuned intermediate 
state) via the absorption of a photon lui and the emission of photon 1^2 then 
Aem = (^1 — W2 — 1^21, and i^{v,t) is equal to fiji^i/A where ftj is the Rabi 
frequency associated with the ujj radiation field, and A is the detuning from 
the intermediate state. Our main interest here is in cw coupling, and so Q{r, t) 
will be independent of time. The key difference between RF coupling and 
Raman coupling is that for RF coupling, the momentum kem is negligible and 
can be ignored. Typically RF coupling (which may extend into the microwave 
region) is single photon, but Raman coupling necessarily involves two photons. 



We note that in Eqs.(44) and ( [45D we have neglected the light shifts that would 
occur for the two photon case, under the assumption that they would be small 
compared to the potentials Vj or the coUisional mean fields near the trap. 



The coupler equations (^^ and (45) can be solved numerically, and we 
illustrate in Fig. |^ a case for an RF coupler in two dimensions, where the 
system begins in an eigenstate of the trap Vi, the RF field is plane wave, and 
V2 = 0. The first two frames show I'f/'ip and |V'2p soon after the coupler has 
been turned on, and we note that the trapped state is essentially unchanged 
from its initial state, while the output state (more or less) copies the initial 
state. After some time, the output beam streams downwards under gravity in 
a well directecLuniform beam that resembles the experimental results of the 
Munich groupE3. 
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Figure 7. Population densities in states |1) and |2) under cw output coupling in a gravita- 
tional field, (a) state |1) at t = 0.424; (b) state |2) at t = 0.424; (c) state |2) at t = 2.544. 
Densities are plotted logarithmically. RF field is plane wave and parameters are Q = 0.2, 
Aem = 40, G = 5, C = 200, V2 = 0. 



7.2 Analytic treatments 

Many of the features of the output coupler behaviour that are found in sim- 
ulations of Eqs.(^4|) and (^) can be understood in terms of analytic approx- 
imations. The key parameters are the strength of the coupling field (il), the 
net energy match with the outpU|t,states, and the bandwidth of the coupling. 
When is large (more preciselyc3 Q'^ + A'^ ^ -I- 2|A|), where fj, is the 
chemical potential of the initial trapped state), it dominates the remaining 
terms in Eqs.(|4^) and (^5|), and the resulting behaviour is synchronous Rabi 
cycling where each part of the condensate oscillates between being fully in 
state |1) and fully in state |2), at frequency fi. The atoms spend such short 
periods of time in the untrapped state, that they are unable to escape. On 



the other hand when is small compared to /i, other terms in Eqs.(44) and 
( p5| ) become more significant, and the difference between the trap potentials 
becomes important. We can understand this by transforming Eqs.(^) and 
( p5| ) into a rotating frame using "02 — '02 exp(iAemt), neglecting kinetic en- 
ergy (because t is small and the atoms do not have time to move), and then 
recognising that the equations simply become Rabi cycling with a spatially 
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dependent effective detuning 

Aeff (r) = Ae,„ - [Vi ~ V2) . (47) 

The mean field contribution to Aes cancels (for w — \) because the atom feels 
the same mean field in either state |1) or |2). Now, because the spectral width 
n of the Rabi cycling is small, significant transfer from |1) to |2) will only 
occur in a small spatial region near the value of r for which Aeff(r) = 0. The 
resulting spatially selective output couphng can be seen in Fig. |7| (b) and (c), 
where the preferred region of coupling is a narrow crescent near y = —10. The 
location of this region can be calculated by noting that the magnetic traps (i.e. 
Vi and V2) are centered on a; = 0, y = 0, and thus the condition Acff(r) = 
for a given choice of Agm is equivalent to selecting an equipotential line for the 
potential difference Vi — V2, which is a circle centered on x = 0,y = 0. Due 
to gravity, the equilibrium centre of mass position of condensate tf^i drops 
below the centre of the magnetic trap to a position where the spring force 
is equal to the gravitational force, as shown in Fig. 0(a). Thus the spatially 
selected region of coupling, where the trap equipotential intersects the trapped 
condensate, takes a crescent shape. Of course the coupling time Tc must be 
long enough that transient broadening does not obscure the spatial selectivity. 



7.3 Weak coupling regime 



The weak coupling regime is of most interest for cw atom lasers, and pertur- 
bative solutions brye been derived which describe key aspects of the output 
field behaviouiEjCZl. Following those treatments we make the assumption of 
an undepleted trapped state and a very dilute output state, which allows us 
to replace ■i/'i in Eq.(^5|) by its initial value il)i{t = 0) and to neglect the 
meanfield interaction of the output state with itself. Taking a Thomas-Fermi 
approximation for V'i(^ = 0), Eq.(E^) becomes a linear Schrodinger equation 



dt 



= -iH, 



cip2 — 7: — WT^i? exp(— l/lt) 



(48) 



where Hoc = — + Voc ■ The potential for the outcoupled wave is Voc = 
V2 + w (/I — Vi ) + Gy and the final term in Eq. (^8|) is a source term in which the 
Thomas-Fermi wavefunction is ^fTF — [l^ — ViY^'^. The eigenfunctions 4>e{i^) 
of Hoc provide a convenient basis for solving Eq.(48). Taking for simplicity 
the one dimensional case (and neglecting gravity) we obtain the solution 



dE I e^s-^^)* 



sin Bt 



B 



(49) 
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where B = [E — + Aem)]/2 and 



A{E) = - J dx(j)*E{x)n{x)exp{ike,nx)'i/TF , (50) 

is the matrix element of the EM coupling amplitude between the output eigen- 
function and the Thomas- Fermi trap state. Eq.(^) provides useful insight 
into the properties of the output field. The term in square brackets is a fa- 
miliar factor from perturbation theory, and at small times the term in braces 
{} — > t so that 'ip2{x,t) « it[^l{x)/2] exp{ikemx)'^TFix). Thus at early times 
the weak coupler produces an output field that is a copy of the initial trapped 
state, spatially masked by the EM field amplitude. At long times the term in 
square brackets [] — > 2'k5{E — fi — Aem) so that the output state becomes 
proportional to the eigenfunction 4>e(x) 

ip2ix, t) « i2Tiex]i{~iE)tK{E)(l)E{x) (51) 

where E conserves energy according to 

= /X + Ae™ , (52) 

and the function A(_E) plavs the role of a spectral filter. The low energy cutoff 
of A(£') can determinecCa by consideration of the classical turning points of 
the wavefunction (j)E{x), which shows that A{E) = for < 0. Thus (for the 
case V2 = 0), the lowest possible output energy for an atom is £' = 0, which 
gives a corresponding lower limit for the coupling field detuning Aem = — A*) 
[see Eq.(p^)] at which value the EM field absorbs all the energy fj, of the 
trapped condensate atom. It is more difhcult to find the upper cut-off of A(i?), 
but the spatial selectivity criterion given by Eq.(|47|) suggests that output 
coupling will cease when Aem exceeds the largest value of Vi — V2 in the 
condensate. Setting V2 — 0, neglecting gravity, and using the Thomas-Fermi 
radius for the edge of the condensate, we can then write the limits for which 
Aem produces effective coupling as 

- H < Aem < ^J■ ■ (53) 

The bandwidth of the output coupling Aujb, which is determined primarily 
by A{E), determines the time scale for the output field to reach its steady 
state (^ [AtJs]^^)- For the case of a Raman coupler (for which e'xjp{ikemx) 
must be retained in Eq.(|50|)) or a gravitational field, Acob is appreciably larger 
than for the case of a simple RF coupler in free space. When fl ^ 
thereby Rabi cycling is suppressed, the faster untrapped atoms escapd£3. 
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7.-^ Finite Temperature Output coupling 

The Oxford group have extended the theory of output couphng to include the 
description of finite temperature condensatesL3E3 (see the article by Burnett 
in this volume). Their treatment, which uses operator fields with a number 
conserving approach is sophisticated, but has an underlying structure that 
is very similar to that presented in the previous section. The atom has two 
internal states \t) and |/) (trapped and free) and the atom field operator is 
divided into 4't{^,t) for the trapped state and ■(/'/(r, for the output field, 
(corresponding to and 'ip2)- The same initial approximations are made in 
assuming the output coupling is weak so that the trapped field is undepleted 
and the output field is very dilute. Inelastic collisions are neglected, and an 
EM coupling field with amplitude given by Eq.(^) models either RF or Ra- 
man coupling. This allows an equation to be written for dijjf{r)/dt, which is 
completely analogous to Eq. (^sj) . The output field moves in a potential com- 
prised of the trap potential for state | /) and the meanfield of the (unmodified) 
initial trapped state, while the self field oi ipf is neglected. The source term 
involves the equilibrium field operator for the trapped condensate, '4>t{Y,t = 0) 
but it is important to note that this is not a meanfield, but includes a ther- 
mal component, in which the quasiparticle levels are assumed to be occupied 
according to a Bose-Einstein distribution of temperature T and chemical po- 
tential /i. The operator V'/ is expanded in modes (/'k(r) of energy that are 
eigenfunctions of a Hamiltonian that is equivalent to -ffoci and leads to the 
solution for the free field 

Vi/(r,t) « ~imv,t)a^{{)) 

-ao(0)^5] [vl/f (r,t)d,(0) + vl/f (r,i)a](0)] , (54) 



where So is the annihilation operator for the condensate mode, TVq is the 
condensate occupation operator ajao, and ctj is the annihilation operator for 
the j*'' quasiparticle level, of energy fi + Ej. The mode functions in Eq.(^) 
are 

^fjir, t)^J2 Mr)Di.,it)Ak,e-'^-* (55) 

k 

for 7? = 0, j+, J- . The quantity D^^it) = t{e-^iK.^' - 1)/(£;1, - E^) 
is an energy selective function in which 

= ^i + A,„, -I- E^ , (56) 
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with Ej± — iEj , and Eq = 0. The quantity Ak,j is a frequency filter analo- 
gous to A(£') in Eq.(|50|) and is given by 

Ak,, = ^ J d^ripl{r)n{x)exp{ike,nx)£,,^{r) (57) 

where for rj = 0, j+, j—, represents V'o, uj, v*, the trapped condensate 
wavefunction, and the usual quasiparticle amplitudes, respectively. We can 
now interpret the three terms in Eq.(p4[). The first term is the condensate 
output and corresponds to the weak coupling mean field calculation given by 
Eq . (^qI) ; the second term {'^■'~^) represents an atom ejected from a trapped 
quasiparticle mode (stimulated quantum evaporation of the thermal excita- 
tions) and the third term (^'■'~) represents an atom ejected from the con- 
densate and creating a trapped quasiparticle (pair-breaking) . Each of these 
processes obeys energy conservation given by Eq.(56), which for the conden- 



sate term {rj = 0) is identical to Eq.(52). For the stimulated evaporation of a 



quasiparticle (77 = j+) we see that the output atom has energy equal to that 
of the initial trapped quasiparticle plus the energy provided by the EM field. 
For the pair breaking term, an atom comes out of the condensate, but leaves 
behind an extra quasiparticle, so that the output energy is /i — Ej plus the 
energy provided by the EM field. 

The density of output atoms is given by 

nout{r,t) ^ {^{r,t)i,f{r,t)) (58) 

and by assuming that all cross correlations between different operators 
ao,aj,a^j vanish and only the populations (ajao) = uq, (ci]dj) = nj, and 
(ajOj) — Uj + 1 remain, we can use Eq.(^J) to write the output density 



nout(r,t) = no|*°(r,t)|2 -1-^ [n,|vl/f (r,t)|2 + (n, + 1)1*^ (r,t)|2j , (59) 

j 



where the three separate components of the output field are easily identified. 
The time evolution of the output field is governed by the output bandwidth 
Aujb as discussed in the previous section, and first discussed extensively by 
Choi et alE3. In the weak coupling limit, fl <C Aw^, we find once again that 
for short times the output state copies the trapped state ^,,(r) 

^J{r,t)^t^o [n{r)/2]Cn{r)t. (60) 

For long times, eu£Tgy conservation selects only output states for which E'k = 
Choi et altB have given a calculation for RF output coupling in one 
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Figure 8. Population output rate from a finite temperature condensate witfi an RF cou- 
pler. Solid line, condensate component. Dashed line, stimulated quantum evaporation 
component. Dotted line, pair breaking component. Parameters: T = 0.5Tc (150?iQj/fc), 
^ = 2.3ftcj/fe , Nt = 2000 atoms. 



dimension, from a condensate at temperature T = 0.5Tc. They calculate the 
rate dNout/dt in steady state, where 



and the result is presented in Fig. |8|, where the three components of the 
output are plotted (as different line types) as a function of the EM detuning 
Aem- A key feature of the output coupling that is readily apparent in this 
figure is how the detuning A^m controls the character of the output field. For 
Aem < —fJ-, only the thermal quasiparticles are ejected from the trap, as we 
expect from energy conservation (see Eg. (|56|) ) : only the quasiparticles have 
sufficient energy to absorb this amount of energy from the radiation field. For 
Aern > —fJ-i up to Agm ~ /i, fuUy Coherent condensate dominates the output 
(see also Eq.(|53|)). Finally for Aem >~ — both pair breaking (where the 
energy of the initial condensate atom plus the photon can be left in the trap) 
and thermal evaporation can occur. 

As a final comment, we note that the processes of stimulated ther- 
mal evaporation_,and pair breaking reduce the coherence of the output 
field. Choi et alE3 have calculated the coherence functions g^-'-'^ (xi, X2,t) and 
g^^"^ {xi, xi,t) and shown that even in the condensate dominated regime, the 
first order coherence is degraded, and g*-^^ exhibits symptoms of atom bunch- 




(61) 



mg. 
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